Quantum transport equations for low-dimensional multiband electronic systems. I 
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A systematic method of calculating the dynamical conductivity tensor in a general multiband 
electronic model with strong boson-mediated electron-electron interactions is described. The the- 
ory is based on the exact semiclassical expression for the coupling between valence electrons and 
electromagnetic fields and on the self-consistent Bethe-Salpeter equations for the electron-hole prop- 
agators. The general diagrammatic perturbation expressions for the intraband and interband single- 
particle conductivity are determined. The relations between the intraband Bethe-Salpeter equation, 
the quantum transport equation and the ordinary transport equation are briefly discussed within 
the memory-function approximation. The effects of the Lorentz dipole-dipole interactions on the 
dynamical conductivity of low-dimensional sp a models are described in the same approximation. 
Such formalism proves useful in studies of different (pseudo)gapped states of quasi-one-dimensional 
systems with the metal-to-insulator phase transitions and can be easily extended to underdoped 
two-dimensional high-T c superconductors. 
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I. INTRODUCTION 

The ab initio band structure calculations of strongly- 
interacting low-dimensional (quasi-one-dimensional 
(Q1D) or two-dimensional (2D)) systems (high-T c 
superconductor YE^CusOr-a; as well as charge- 
density-wave (CDW) systems K0.3M0O3 and BaVS3 
being examples 1 - - — ) usually reveal several bands in the 
vicinity of the Fermi level that are, in the leading 
approximation, decoupled from the rest of the band 
structure. In the tight-binding approximation, the bare 
dispersions of such a multiband model are well under- 
stood in terms of the strong hybridization among several 
properly chosen orbitals per unit cell, while the hy- 
bridization with the other electronic levels is neglected; 4 
The bare single-particle tight-binding Hamiltonian is 
thus taken to include commensurate/incommensurate 
CDW or SDW (spin-density wave) average fields and 
several first-neighbour bond energies. In the two-particle 
Hamiltonian one usually retains the long-range Coulomb 
interactions, the largest local (Hubbard) interactions to- 
gether with the most important short-range interactions. 
Associated with different terms in the single-particle 
Hamiltonian are fluctuations in the charge and spin 
densities on different orbitals in the unit cell and on 
the corresponding bonds. All these ingredients make 
strongly-interacting multiband electronic problems 
extremely complicated. 

To explain most of experimental results on typical 
single-band systems, it usually suffices to describe prop- 
erly the coupling between q ss monopole charge fluc- 
tuations and external electromagnetic (EM) fields and 
combine this coupling Hamiltonian with the ordinary 
transport equations. 5,6 However, in multiband electronic 
systems, even in the cases when only one band inter- 
sects the Fermi level, the analysis is considerably more 
complicated. Not only the interband physics but also 
the intraband physics depends drastically on both single- 



particle and two-particle parameters of the Hamiltonian 
and are affected differently by various types of fluctu- 
ations. Therefore, to understand the results of trans- 
port measurements 7 -^ and reflectivity measurements*^— 
of multiband electronic systems in detail, we describe the 
q«0 multipolar fluctuations in the dynamical conduc- 
tivity tensor in a way consistent with the results of experi- 
mental methods which probe charge and spin fluctuations 
which presumably occur at finite wave vector o 12 ' 13 . The 
q«0 charge-charge correlation function is thus just one 
of several equally important correlation functions defined 
in terms of symmetrized intraband and interband charge 
and spin fluctuations. Common to all these correlation 
functions is a complicated structure of the intraband and 
interband electron-hole propagators. Therefore, to un- 
derstand the contribution of different types of fluctua- 
tions to these correlation functions, in particular to the 
dynamical conductivity, we must first find the solution to 
quite general Bethe-Salpeter equations for the electron- 
hole propagators and then combine the result with the 
very definition of the response function under consider- 
ation. This question, which is of great importance in 
studies of strongly-interacting low-dimensional systems, 
is in the focus of the present analysis. 

In this article, we thus present a general response 
theory which uses both an essentially exact semiclas- 
sical description of the coupling between valence elec- 
trons and external EM field s 14 ! 15 and a quite gen- 
eral description of the intraband and interband charge 
fluctuation s 16 ! 17 . The result is the current-dipole Kubo 
formula for the dynamical conductivity^ with the struc- 
ture of the intraband electron-hole propagators deter- 
mined by means of the quantum transport equation (i.e., 
the intraband Bethe-Salpeter equation). To emphasize 
the connection between the present theory and the or- 
dinary transport theory,— &22r2L we determine the ex- 
plicit form of the intraband conductivity in the memory- 
function approximation^ - — We also find the structure 
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of the interband conductivity in both cases, with real 
and with purely imaginary interband dipole vertices, to 
show in which way the Lorentz local held corrections^— 
affect the dynamical conductivity of interacting low- 
dimensional systems. 

The article is organized as follows. To make the read- 
ing of the article easier, we give in section 2 an overview 
of the results. Then in section 3 the notation is in- 
troduced to describe electrons and boson modes in a 
general multiband electronic model with boson-mediated 
electron-electron interactions. A complementary multi- 
band electronic model with nonretarded electron-electron 
interactions will be studied in a separate article 2 ^ which 
is focused on the commensurability effects in the single- 
particle conductivity tensor. 

Starting with the quantum electrodynamical descrip- 
tion of the interaction between valence electrons and ex- 
ternal EM fields, we derive the exact semiclassical form of 
the electron-EM field coupling Hamiltonian in section 4. 
In section 5 we define the dynamical conductivity ten- 
sor and briefly discuss the relaxation-time approxima- 
tion. Section 6 gives the structure of the self-consistent 
Bethe-Salpeter equations for the dipole vertex functions 
and for the related auxiliary electron-hole propagators 
(to be defined latter) in a general multiband case. We 
show the equivalence of the Bethe-Salpeter equation for 
the auxiliary intraband electron-hole propagators and the 
quantum transport equation. The generalized Drude for- 
mula is derived from the quantum transport equation by 
using the memory-function approximation. The bare in- 
terband conductivity is considered in section 7 and the 
related local field effects in section 8. Section 9 contains 
the concluding remarks. 



II. OVERVIEW OF RESULTS 

Before discussing the details of the dynamical con- 
ductivity model proposed in this article, we present an 
overview of the results. In this way principal physical 
messages of the present work are separated from the tech- 
nical details. 

We are interested here in a general low-dimensional 
model with multiple electronic bands in the presence 
of boson-mediated electron-electron interactions. The 
single-particle part of the total Hamiltonian is assumed 
to be exactly solvable. The Bloch energies £z,(k) in the 
band labeled by the band index L and all relevant bare 
vertex functions (the dipole vertices P^ L (k, k'), for ex- 
ample) are thus taken as known functions. To study the 
influence of the boson-mediated electron-electron inter- 
actions on single-electron propagators ^(k, iw n ) and, in 
addition, on different electron-hole propagators encoun- 
tered in the response functions of interest, we use high- 
order diagrammatic perturbation theory of the metallic 
state, as usual in investigations of low-dimensional elec- 
tronic systems i 29 ' 30 

The electron-hole propagators are associated with the 
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FIG. 1: The Dyson equation ([TJ for electron Green's func- 
tions. L is the electron band index and v is the boson branch 
index. The black circle and the big gray circle represent, re- 
spectively, the bare and the renormalized electron-boson ver- 
tex. Here and herafter, the tiny and bold solid (dashed) lines 
are the bare and renormalized electron (boson) propagators. 

dispersive, undamped bosons which can either corre- 
spond to the external degrees of freedom coupled to elec- 
trons/holes or the internal degrees of freedom built of 
the electron-hole excitations at finite values of q. In sec- 
tions 6-8, the emphasis is on the acoustic and Raman- 
active and infrared-active optical phonons. The scatter- 
ing by phase and amplitude phonon modes of the com- 
mon Peierls CDW model will be studied in detail within 
the microscopic Lee-Rice- Anderson model in the accom- 
panying article 3 -^ (to be referred to as Article II). 

In high-order perturbation theory, the electron Green's 
function Si(k, hd n ) satisfies the Dyson equatio n 29 i 32 ' 33 

[ifiwn - £i(k) hY. L (k,iuj n )]g L ('k,iu> n ) = H, (1) 

where Sj,(k, iui n ) is the self-energy function of the elec- 
tron (sec figure 1). In this case Si(k, iu) n ) is a func- 
tion of £/i(k, iu) n ), resulting in a system of two coupled 
equations. A strong reduction of the quasi-particle pole 
weight in (?z,(k, iw„) at the Fermi level is predicted by dif- 
ferent strictly ID theorie s 30 ' 34 ' 35 and observed in angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ments on Q1D system a 36 ' 37 . In this article we present 
a systematic method of calculating the dynamical con- 
ductivity tensor in interacting low-dimensional systems 
which is consistent with the described treatment of the 
single-electron Green's functions Sz,(k, iuj n ). 

A. Electron coupling to EM fields 

In semiclassical electrodynamics applied to a general 
exactly solvable single-particle problem with multiple 
electronic bands, the electric component of the EM field 
parallel to the cartesian axis a, E a (q), couples to the cur- 
rents of free charges and to the dipole moments of bound 
charges in a way described by the coupling Hamiltonian 

H ^ = -J2E a (q)P a (-q). (2) 

qa 

Here, P Q (— q) is the total dipole density operator. It is 
the sum of the intraband and interband contributions 

£«(-q) = ^^P a L,I (k + ,k)4 +q/ikff , (3) 

LL' kcr 
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with L (k + , k) representing the corresponding intra- 
band [L = V) and interband (L ^ L') dipole vertices, 
and k + = k + q. It is not hard to verify that these ver- 
tices are directly related to the vertices in the current and 
charge density operators, J„ L (k + ,k) and q L L (k + , k) . 
The relation is given by the expression 

(q = q a e a in the longitudinal case). The imaginary cur- 
rents and the imaginary dipole moments do not have a 
straightforward physical meaning, but they make the no- 
tation used throughout this article as simple as possible. 
For example, they can be used to express the gauge in- 
variance of the expression @ simply by rewritting iJ ext 
in terms of q L L (k + , k) and L (k + , k) and using the 
standard relations between E a (q, lu) and the scalar po- 
tential V tot (q,uj) and the vector potential A a (q,ui). 



B. Dynamical conductivity tensor 

In order to examine how electrons respond to ap- 
plied EM fields, we introduce the dynamical conductiv- 
ity tensor a a p(q,uj)r^ which is the response function re- 
lating the induced total current J Q (q, ui) to £73 (q, w); 
J Q (q, ui) = ""^(q, w)Ep(q, lu). We use here the fi- 
nite temperature formalism. The dynamical conductivity 
cr a| g(q,£j) is in this case given by analytical continuation 
(iv n — ¥ lu + ir/) of cr Q( 3(q, i^„), which is the Matsubara 
Fourier transform of <r aj g(q, t). The compact form of the 
semiclassical coupling Hamiltonian @ allows us to ex- 
press a a j3 (q, r) in the form of the current-dipole Kubo 
formula^ 

Vapfar) = ■^ 7 (T r [J a (q,r)P (3 (-q,0)]) irrcd . (5) 

Here, J a (q) is the total current density operator defined 
in a way analogous to ([3]) (see (I^Bl ). In order to extract 
the most important physical consequences that follow 
from this expression, one usually uses the relations (U) 
to rewrite a a p (q, r) either in terms of the charge-charge 
or the current-current Kubo formulae^ In contrast to 
these approaches, we use here JS]), because this expres- 
sion will be shown below to provide the direct link of the 
present analysis to the common transport theory^ and 
to the textbook approaches to the dynamical long-range 
Coulomb screenin g 32 i 33 i2&2£. 

To complete the response theory, we must take into 
account the relaxation processes describing the way in 
which electrons are scattered by static disorder, by other 
electons or by different types of boson modes. Since 
our primarily interest is in the low-dimensional sys- 
tems where the quasi-particle pole weight in the single- 
electron Green's function is strongly reduced by boson- 
mediated electron-electron interactions, we must calcu- 
late a a (j (q, t) beyond the common weak coupling the- 
ory (corresponding to second-order perturbation theory 
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FIG. 2: The Bethe-Salpeter expression for the screened dy- 
namical conductivity tensor a a p(f\, if n ) in a general multi- 
band case. The dynamical conductivity a a p{<\, iv n ), equation 
([5]), is the RPA-irreducible part of cf a p(q, i«/„).— 



for the single-electron self-energy function). This means 
that in order to work out a a/ 3 (q, r) we must use the 
general Bethe-Salpeter expression for the screened con- 
ductivity tensor illustrated in figure 2 and the definition 
Y,p°o t p{<i.,u)Ep(q L ,u) = Y,i3°ai3(<l,oj)E p(ci,uj), where 
Eo(q, cv) is the external electric fielcU^. The bold solid 
lines in these two diagrams are the single-electron Green's 
functions which satisfy the Dyson equation ([T]). The in- 
traband and interband contributions to a a p (q, if„) = 

Sll' a aB (Qji^n) can be represented by the same ex- 
pression, 



^|'(q>n) = ^E^2 E ^ Li '(k)SL(k,k,„) 



xg L ,(k + ,iuj n+ )(-)Tp L (k,k + ,iw„,kj„ + ), (6) 



equally well for L = L 1 as for L 7^ L' . Here the 
Tp L (k, k + , iw„, iLu n+ ) are the renormalized dipole vertex 
functions which satisfy the self-consistent Bethe-Salpeter 
equations for the renormalized dipole vertex, and \Lu n+ = 
iuj n + \v n . In order to emphasize the connection between 
the present approach and the common transport theory, 
as well as to simplify the notation, we introduce the aux- 
iliary electron- hole propagator (k, k + , iw„, iuj n +) by 
comparing the relation 
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Jf L (k)$^ L (k,k+>„>„+) 

(7) 

with ©. Therefore, our task is to solve the self-consistent 
equations for $a (k, k + , iu„ , ilu 71+ ) rather than that for 
the renormalized dipole vertices. 

We limit our discussion to the low-dimensional systems 
where the scattering by static disorder and the scattering 
processes in which electrons change the band can both be 
neglected. These processes have negligible influence on 
the properties studied here. However, the self-consistent 
equations still have a very complex structure in particular 
as a consequence of the following features: the nesting 
properties of the Fermi surface, the resonant nature of the 
scattering by the (quasi)static CDW or SDW potentials 
and the local-field effects in the interband conductivity 
channel. 
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C. Results 

The results of the present analysis are the following: (i) 
we show the formal solution to the Bethe-Salpeter equa- 
tions for <&^l (k, k + , icj„, iuj n+ ) and to the (intraband) 
quantum transport equation, (ii) we compare the intra- 
band conductivity to that obtained by using the memory- 
function approximation, (Hi) we calculate the bare inter- 
band conductivity and (iv) we rederive the well-known 
Lorentz-Lorenz form of the optical conductivity in the 
presence of scattering by boson modes. 

(i) Although it is tempting to completely neglect the 
electron-EM field vertex corrections in a a p(q, this is 
not justified in the intraband channel, at least in the cases 
with low-dimensional conductivity studied here. The 
method by which the intraband conductivity tensor is 
calculated here is very general and is based on a sys- 
tematic analysis of the self-consistent equation for the 
auxiliary electron- hole propagator ^(k, k + , ioj n , \uj n+ ) 
in which the single-electron self-energy contributions and 
the related vertex corrections are treated on an equal 
footing. The equation is of the form 

[ii/ n + E(k)/S-e(]q.)/fi 

+n(k, k +! iw n , ioJ n+ )] $ a (k, k + , iw„, ioj n+ ) 



L'j k'+q L' k+q L\ k'+q L' k+q L\ k'+q L' k+q 



— [g(k,iu n ) - g(k + ,iu n+ )]p a (-q) 



(8) 



with p (— q) = ie/q a and II(k, k +) iu) n , iw„+) = 
n r (k,k+,iu;„,ia; n+ ) + iIP(k, k+, iuj n , iw n+ ). For sim- 
plicity we assume that only one band intersects the 
Fermi level and omit explicit reference to the conduc- 
tion band index. II(k, k + , iui n , iuj n+ ) is the electron-hole 
self-energy, the structure of which is derived here for the 
case in which conduction electrons are scattered only by 
various boson modes. These boson modes are assumed 
to dissipate momentum by their own means (impurities, 
Umklapps, etc.). Pinned phasons of the commensurate 
CDW problem, for example, enter in the present theory 
in the same way as the infrared-active optical phonons 
from section 8.2. This question is discussed in more de- 
tails in Article II. 

The electron-hole self-energy is found to be the differ- 
ence between two auxiliary single-electron self-energies 
E(k, Lv n ), and the contributions of the q rj for- 
ward scattering processes are found to drop out of both 
S(k, iui n ) and II(k, k + , iuj ni hu n +), in full agreement with 
the charge continuity equation. 

Equation © is valid under quite general conditions, 
and is therefore applicable to a rich variety of problems, 
including purely electronic models with nonretarded local 
and short-range interactions. For example, it is not hard 
to verify that the normal scattering processes do not lead 
to the resistivity in this case, as required by the general 
principle of microscopic reversibility^ This question will 
be discussed in more details in Ref. l28l 

(ii) In the ordinary 3D metallic regime, equation ((5J) 
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FIG. 3: The Bethe-Salpeter equations for electron-hole prop- 
agators in a general multiband electronic model with weak 
nonretarded and/or boson-mediated electron-electron inter- 
actions. The vertex renormalizations are assumed to be neg- 
ligible in this case. 



leads to the ordinary transport equation, with 5 ^ 

.intra, . , ™Z gO ( ~ gWggW) (q) 

being the resulting intraband electrical conductivity, and 
q = q a e a again. Here, m is the unity of mass, com- 
monly taken as the free electron mass, and Fi Q (k) oc 
IT(k, k+, \w n , iui n +) is the intraband relaxation rate. 
In the Drude uj rs limit, where uj 3> q a v a (\i), 
the q a dependent term in the denominator of © can 
be ignored, and we obtain the intraband conductiv- 
ity (fTUl) . For uj <C q a v a (\i) and q a —> 0, on the 
other hand, the result is the well-known Thomas-Fermi 
expressio n 33 ' 38 for c™* ra (q, uj) and for the dielectric sus- 
ceptibility X intra (q,u;) = (^/iwjo^'fq, w). 

(Hi) The general expression for the interband single- 
particle conductivity is derived and compared to the ap- 
proximate expression which is frequently encountered in 
the literature i^ 8 -^ 

(iv) If the scattering of valence electrons is sufficiently 
weak, the electron-boson vertex renormalizations are of 
no importance and the electron-hole propagators of fig- 
ure 2 satisfy the Bethe-Salpeter equations shown in 
figure 3 with two types of contributions, the ladder 
contributions and the RPA (random phase approxima- 
tion) contributions^ These equations can be completely 
rewritten in terms of the renormalized single-electron 
Green's functions and the dynamically screened electron- 
electron interactions. 

This approximation is used here to illustrate how the 
conductivity tensor is affected by the short-range dipole- 
dipole interactions. The result is the well-known Lorentz- 
Lorenz form of the dynamical conductivity for 3D metal- 
lic sp a models 25 '-- which is believed to be valid in the 
Q1D conductors, as well2£. From the result 
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FIG. 4: The Dyson equation {16} for boson Green's functions. 

we can conclude that the contributions to the conductiv- 
ity tensor of the free charges and of the bound charges are 
decoupled and that the Lorentz dipole-dipole interactions 
affect only the bound (interband) charge contribution to 
the conductivity tensor. 



III. ELECTRON-BOSON COUPLING 

Let us first introduce the quantities characterizing a 
general multiband electronic model with boson-mediated 
electron-electron interactions. The total Hamiltonian is 
the sum of two contributions, H = Hq + H'. Here, 

H Q = < + = EL(k)4 ka C Lk . 

Lkcr 

+ E 2M: \- p U'P»i' + ( M ^q') 2 <4'<w] (12) 

l/q' 

is the bare electron-boson Hamiltonian, and 

ff' = EE <&r L (k + q',k)4, k+qV c LkCT (13) 
z^q' LL'kir 

is the scattering Hamiltonian which describes the scatter- 
ing of valence electrons by each of the boson modes. As 
already mentioned, the bosonic modes can be external to 
electronic system (e.g. phonons) or built from electron- 
hole fluctuations at finite wave vectors. In (fl~2|) . £i(k) 
is the bare electron dispersion, oj^q' is the bare boson 
frequency and M v is the corresponding "mass" parame- 
ter. Finally, $J L (k + q', k) is the generalized scattering 
potential. For convenience, it is assumed to be a prod- 
uct of the coupling constant g u /\^N, the dimensionless 
electron-boson vertex L (k + q', k) and the boson field 

<^' L (k + q',k) = q L J L {k + q\ k) W . (14) 

In the intraband scattering approximation, in which the 
electron does not change the band when it is scattered 
by bosons, the scattering potential $^ L (k + q', k) is di- 
agonal in the band index. 



A. Skeleton diagrams 

In the finite temperature formalism, the effects of H' 
on the single-electron and boson properties are described 
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FIG. 5: The skeleton series for the single-electron self-energy. 
The self-energy is shown in terms of both the renormal- 
ized three-point vertex function L (k+, k, kj n +, ii^n) (the 
big gray circle) and the completely irreducible four-point ver- 
tex function U„ L,LL (k + q', k, q', \u n +, iui n , iv n ) (the white 
rectangle). 

in terms of two exact Matsubara Green's function o 29 ' 32 ' 33 

fo(k,r)=-(T T [ CLkCT (r)4 kCT (0)]) ! 
^(q',r) = -(T T [ W (r)gt q ,( )]). (15) 

Their Fourier transforms satisfy two Dyson equations, 
equation ([1]) and 

{M v [{iv n f-u 2 v ^] -fiE, / (q',ii/ n )}X> v (q',ii/ n ) = H (16) 

(see figure 4). It is convenient to adopt an abbreviated 
notation by writting 

^^(k+q'^l^qV) EE ftjf L (k+q',k,r). (17) 

This expression can also be shown in the form T„ L (k + 
q',k, t) = J"^' L > LL '(k + q',k, q',r), where 

^'^^(k + q'.k'^r) 

= -<T r [^' L {k + q', k, r)^ lL[ (k', k' + q', 0)] ),(18) 

is by definition the force-force correlation functio n 33 i 41 
(FFCF) associated with the potential <fr^' L (k + q', k). 

In high-order perturbation theory the self-energy func- 
tion of the electron, Ei(k, kj„), is usually represented by 
the sum of skeleton diagrams shown in figure 5, leading 
to a simple analytical expression 

hY, L (k,iu n ) rj - ■^E^ i '( k + q '' ia;n +) 

x £L'L,LL' {k + q / k q ' ; iWn+; iWn ). (19) 

Here, T„ L ' LL (k + q', k, q', iw n+ , iuj n ) is the renormal- 
ized FFCF obtained by replacing one bare electron- 
boson vertex in F^' L ' LL '{k + q', k, q', iv n ), q(;' L (k + ,k), 
by the renormalized vertex L (k+, k, icj„ + , iu) n ). Equa- 
tion (|19|) is nothing but the three-point vertex representa- 
tion of E^(k, iw n ). In the same representation the boson 
self-energy takes the form 




x q h v h ' (k, k+ )Q L (k, iu n )Q L < (k+ , ioj n+ ) . (20) 
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We have just reemphasized that the three-point vertex 
representation is very useful in describing the closed sys- 
tem of four equations characterizing strongly-interacting 
electron-boson systems in the absence of external EM 
fields (equations {TJ), tfl6]>. (fT9)l and ([20])V However, the 
conductivity tensor of such an electron-boson model has 
a compact form if the structure of the Bethe-Salpeter 
equations of section 6 is also simple, at least from the 
formal point of view. For this purpose, it is much 
more convenient to use the four-point vertex representa- 
tion. This representation is based upon the knowledge 
of the irreducible four-point interaction U% L ' LL (k + 
q', k, q', iw n+ , iw n , iv n ). The equivalence of the two rep- 
resentations in the skeleton approximation is illustrated 
in figure 5, where U% L ' LL (k + q', k, q', iw„+, iw„, w n ) is 
represented by the white rectangle. The expression for 
S^(k, iu) n ) is now of the form 

hT, L {k,iuj n ) w - ^2 ^pE^ L '( k + q '' iWn +) 



fr 1 , P h z 



xU. 



L'L,LL' 



(k + q', k, q', iu n+ , iw n , iv n ). (21) 



A full discussion of the Bethe-Salpeter equations is post- 
poned to section 6. 



IV. ELECTRON COUPLING TO EM FIELDS 

In quantum electrodynamics of tightly bound valence 
electrons in solids, whether metallic or insulating, the 
coupling between electrons and external EM fields is 
given by the minimal gauge invariant substitution. The 
expansion to the second order in the vector potential 
A a (q,u>) gives the coupling Hamiltonian H° xt ss Hf xt + 
#| xt , withil^^ 



flr , = -i^^(q)J a (-q), 



qa 

^ oxt = ^ E ^(q-q')^(q')7"(-q;2). (22) 



qq'a^ 



Here, 



4(q) = EE J "'( k > k +) c L CT CL'k +qCT , 



LL' kcr 

7"(q; 2) = EE ^ k > k+; 2)c[ kcr c Lk+qff (23) 

L kcr 

are, respectively, the total current density operator and 
the bare diamagnetic density operator. 

To obtain the semiclassical expression for the coupling 
Hamiltonian, equation ([2]), we have to take into account 
the fact that both the interband and intraband current- 
current contributions to the diamagnetic current must 
vanish in the normal metallic or insulating state. The 



cancellation of the interband contribution is a conse- 
quence of the effective mass theorem; 17 i 24 ' 43 i 44 



7"(k) = 7 "(k;2)+ 

_ TO 9 2 £i(k) 

h 2 dk a dk/3 ' 
leading to H cxt m #f xt 

r — ' 



(2m/e 2 )J^ (k)J^(k) 
e L (k)-e L /(k) 



(24) 



H£ xt , withii 



£j Q "(-q) + ^ 

L L^L' 



i^^ Li '(-q) 



rrcxt 
H 2 



2mc 2 



£ A Q (q-q')^(q')7^(-q) 



(25) 



Lqq'ct/3 



Here, P^ L (q), L ^ V , is the interband part in dipole 
density operator ((3J, and 7^3 (— q) is the diamagnetic 
density operator. The reciprocal effective mass tensor 
7^(k) is the vertex function in 7^g(— q)- 

The intraband current-current contribution to the dia- 
magnetic current cancels H^, confirming the assertion 
made in section 2 that 



E 

qO: 



E a (q)P a (-q). 



We emphasize the simple structure of this coupling 
Hamiltonian. The most direct way to obtain this ex- 
pression for H cxt is to consider the coupling of the scalar 
potential V tot (q, uj) to the charge density operator^ 



A. Density operators 

Let us derive now the relations Q among three types 
of intraband and interband vertex functions underlying 
the gauge invariance of the present response theory. To 
do this, we use very general microscopic operator equa- 
tions 

d 

hq a J a (q) = i^/3(q) = [/5(q),#o'L 

ihJ a (q) = i^A*(q) = [P a (q),H$] , (26) 

where is the electronic part in the bare Hamiltonian 
(TT2"1) . Furthermore, 



P(q) = E E e 9 LL ( k ' k +)cL CT C i'k+qr 



(27) 



LL' kcr 



is the total charge density operator, with eq LL (k, k+) 
being the charge vertex function, while J a (q) and P a (q) 
are given by (|23|) and (|3|). The two commutators in (|26|) 
are easily worked out, leading to the required relations 



C (k,k+) 



ifiJ^(k) _ ie „ 



e L '(k+) - e L (k) q a 



-q LL (k,k+) 



The explicit form of each of these vertices is often ob- 
tained by calculating the coefficients in the linear term 
in iJ 0Xt in ([22]) . 

The long-wavelength intraband vertices are of partic- 
ular interest. They are given by model independent ex- 
pressions 

g iL (k,k + )«l, 

J„ LL (k,k + )« j£ £ (k) = et£°(k), 

P^ L (k, k+) » P Q Li (-q) = p a (-q). (28) 

Here, ?;^ ,0 (k) = (l/K)dei,(\s)/dk a is the bare electron 
group velocity and p a {— q) = ie/q a is the bare dipole 
vertex. 

Notice that the thermodynamic averages of the den- 
sity operators in (|26|) are the usual induced macroscopic 
densities. For example, J(q, lu), p(q,uj) and P(q, w) rep- 
resent the common notation for the induced macroscopic 
current, charge and dipole densities. They are related by 
two equations^ 

q- J(q,w) =wp(q,w), 

p(q,w) = -iq-P(cbw). (29) 

The first one is the usual charge continuity equation. 



V. DYNAMICAL CONDUCTIVITY TENSOR 

We now define the total dynamical conductivity tensor 
<r ai g(q, lu) as the response function to the macroscopic 
electric fielcU^ 

E(q,u;) =Eo(q,w) + Ei(q,w); (30) 

that is 

J Q (q,w) = ^o- Q( 3(q,w)P (3 (q,a;). (31) 



Here, J(q, lu) is the induced total current density from 
([29)) . Eo(q, lu) is the external field and Ei(q, lu) is the 
depolarization field (equal to — 47rP(q, lu) for the longitu- 
dinal polarization of the field in the thin slab). The con- 
ductivity tensor is given by analytical continuation to the 
real axis of <J a| g(q, iv n ) (see figure 6), where <j a| g (q, iu n ) 
is the Matsubara Fourier transform of 

Ca/9(q,T) = ^-(T T [j a (q,T)Pp(-q,Q)]) ixmd . 

Therefore, to find the dynamical conductivity tensor of 
an interacting electron-boson system, we must deter- 
mine the structure of the current-dipole Kubo formula 
([S]) which agrees with the definition of the macroscopic 
electric field (|HJ). 
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L' k+q L' k+q LJ k'+q L' k+q 




FIG. 6: The dynamical conductivity cr a /3(q, iv n ) from fig- 
ure 2, shown in the three-point vertex representation (left 
diagram) and in the four-point vertex representation (right 
diagrams). The collection of diagrams in the renormalized 
four-point interaction must be consistent with the definition 
of E(q, lu) in (|30|) ; i.e., the renormalized four-point interaction 
is irreducible in the long-range RPA contribution. 

A. Relaxation-time approximation 

It is instructive first to show the structure of the 
bare dynamical conductivity ([5]) in the simplest multi- 
band model, namely, a two-band model with the lower 
band partially occupied by electrons and the upper band 
empty, in the limiting case where the scattering is suf- 
ficiently weak. As explained in more details in Article 
II, this model is of importance in accounting for the 
properties of the ordered state of some simple CDW 
systemsi 41 i 45 ' 46 

To do this, we first calculate the ideal bare conductivity 
tensor a a/3 (q, \v n ) which is given by the first term in fig- 
ure 6, with Gl(q, iw„) replaced by the bare Green's func- 
tion Ql(q,iuj n ). Summation over icu n can be performed 

easily to give o-°^(q, iv n ) = 6 ai p Y^lv a aa (^' 1L, n) ■ Af- 
ter analytical continuation, we obtain the ideal bare in- 
traband (L = L') and interband (L ^ L 1 ) contributions 
described by 

i v Jt L '(k)pj' L (k)[f L (k) - f L ,(k + )} 

nuj + ih v + e L (k)-e L ,(k + ) • 1 ' 

Weak scattering processes can be inserted into this ex- 
pression by replacing the adiabatic term 77 by the relax- 
ation rate T^ L independent of frequency and wave vec- 
tor. This is the generalization of the common relaxation- 
time approximation. As required, the associated trans- 
port relaxation time = 1/T^ L differs from that found 
in the single-electron propagators. 

The bare intraband contribution is of the Drude form 

intra.O/ , ,\ ^ ie 2 n°^'° /m 

°aa (q> w )~ > (33) 

LU + 11 i a 

with 

n ^ = 7^ m ^)] a (-^) (34) 

representing the effective number of conduction 
electrons ) 17 ' 41 and T^ L = Ti a by assumption. n ac ; 
includes the conduction holes from the lower band and 
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L'j k'+q L' k+q h\ k'+q L' k+q L'jk'+q 
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+ 
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FIG. 7: The general form of the Bethe-Salpeter equations 
for electron-hole propagators. The white rectangle represents 
again the completely irreducible four-point interaction. 



FIG. 8: The Bethe-Salpeter equations (|36[) for auxiliary 
electron-hole propagators. 



the thermally activated electrons from the upper band. 
The bare interband contribution is the sum of two 
off-diagonal terms in (TJ2"1) . with r\ replaced by T^ L : 

< er '°(q^) = E^L^"'(q^)- 

The main disadvantage of this simple procedure of 
calculating er QQ (q, uj) is that it does not apply to low- 
dimensional systems with sizeable interactions. The rea- 
son for that is the following. First, comparison with ex- 
perimental data shows that the relaxation rates T^ L are 
frequency dependent. Second, the causality principle re- 
quires that T^ L (ui) is just the imaginary part of a com- 
plex quantity which will be called here the electron-hole 
self-energy function. Finally, it turns out that in the 
strong coupling limit the electron-hole self-energy is not 
a simple function of frequency ui. Instead, it depends 
quite drastically on two internal electron frequencies too. 

We shall therefore present now a systematic treatment 
of scattering processes in interacting electron-boson sys- 
tems with multiple electronic bands which treats all el- 
ements in the electron-hole self-energy consistently with 
the charge continuity equation (f25)) . with the causality 
principle and with the definition of E(q, uj). The result- 
ing self-consistent equations turn out to be the general- 
ization of the equations used by Vollhardt and Wolflc 21 
to study the scattering by static disorder in a single-band 
case. 



VI. BETHE-SALPETER EQUATIONS 

The term Bethe-Salpeter equation s 40 ' 47 refers to the 
self-consistent equations for the exact electron-hole prop- 
agators <^ll'l' 1 l 1 ( k) k + , k^_, k', iw„, iw„+, iuj m+7 iui m ) 
shown in figure 7. A widely-used simplification in these 
equations is illustrated in figure 3. In this case, we 
have both the ladder contributions and the RPA-like 
contributions to the irreducible four-point interaction, 
but the three-point vertex renormalizations are ne- 
glected both in these equations and in the related Dyson 
equations. The interaction lines and the electron propa- 
gators are renormalized, so that these equations must be 
treated consistently with the Dyson equations, as already 
pointed out in section 2. However, in strictly ID metallic 
systems as well as in strongly-correlated low-dimensional 
systems, it is necessary to use the Bethe-Salpeter 
equations beyond this approximation by taking into 
account the three-point vertex renormalizations as well. 

To simplify the analysis, it is convenient to multiply 



the left-hand side of the diagrams in figure 7 by the cor- 
responding bare vertex function and then sum over k' and 
ibj m . For the electron- hole propagators in the dynamical 
conductivity ([5]), the bare vertex in question is the bare 
dipole vertex function, and the Bethe-Salpeter equations 
from figure 7 turn into the Bethe-Salpeter equations for 
the auxiliary electron-hole propagators 

$^ L '(k,k + ,kj„,iw n+ ) = -^ L (k,iw„) 

x<?L'(k+,iw„ + )(-)r" (k, k + ,kj„,iw n+ ). (35) 

The equations are illustrated in figure 8. 

With little loss of generality, we restrict the analysis 
to the case where the irreducible four-point interaction 

Yj V Uv u 1 (k + ,k',q',iw n+ ,iw m ,n/ m ) can be approxi- 
mated by U L i i' i i i (q' 5 \v m ). In this case, the explicit 
form of the intraband [L = L') and interband (L ^ L') 
equations is the following 

-D~^,(k, k + , iw„, iu; n+ )$" (k, k + , iu n , iuj n+ ) 

= - <?z/(k + ,iw n+ )] 

x{-P Q L ' L (k + ,k)-i]T J2 U L ' L '^ L {^,iv m ) 

^ P L 1 L' 1 q'ii^,, 

X Qa 1 L[ (k + q' , k+ + q' , lUJ m , iu m+ ) J , (36) 

with 

D^,(k,k + ,iw„,ia; n+ ) = iu n + e L {k)/h - £ L >(k + )/h 

+SL(k,iw„) - £ L /(k + ,kj„ + ) (37) 

and iuj rn — iuj n + \v m . The 5i(k, iuj n ) satisfy the Dyson 
equations for electrons ((TJ and U L Ll ' 1 (q',ii/ m ) is a 
simple generalization of the completely irreducible four- 
point interaction [/ L ' L > LL '(q', iv m ) from pi]). 

In the electronic multiband model we are interested 
in here, there are n bands in the vicinity of the Fermi 
level. The bands are built of n orbitals per unit cell. 
Associated with each of these orbitals are the intracell 
charge and spin fluctuations, and these intracell fluctu- 
ations are possible sources of scattering. The monopole 
charge fluctuations underlying the present analysis of the 
dynamical conductivity tensor might be also strongly af- 
fected by the presence of the intracell fluctuations. As 
a result, cr ai g(q, uj) will have considerably more structure 
than that characterizing the ordinary single-band model, 
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even in the case where only one band out of n valence 
bands intersects the Fermi level. 



A. Quantum transport equation 

In the remainder of this section we consider the intra- 
band term in Q36[) from this point of view. We assume 
that there is only one conduction band and again drop 
reference to the band index. We use the intraband scat- 
tering approximation where the scattering processes in 
which the electron changes the band are neglected. As 
discussed in Article II, this approximation is easily seen 
not to be adequate, for example, to study fluctuations 
in the order parameter in the pseudogap state of the or- 
dinary CDW systems. However, it is convenient to il- 
lustrate the main ideas regarding the quantum transport 
equations. Finally, the three-point vertex renormaliza- 
tions in the renormalized interactions are also neglected, 
as usual in the common non-crossing approximation for 
the electron-boson systems4£~— The resulting intraband 
Bethe-Salpeter equation is of the form 

D^ 1 (k, k + , kj„, iw n+ )$ Q (k, k + , iw n , iu) n+ ) 

= [£(k,iu; n ) - g(k + ,iuj n+ )] jp a (-q) 

-4 X! •7 r (q'> iiy m)'I > a(k + q',k + + q',iu> m ,iw m+ )l. 

(38) 

The structure of (|38f is simplified by defining the aux- 
iliary single-electron self-energy function 

^(k, 1 c„)=--^ll-^ r j 

xJ"(q',ii/ m )£(k + q',iu; n + H/ m ) (39) 

and the electron-hole self-energy 

II(k,k + ,iu; n ,kj„ + ) = E(k,iw n ) - E(k + ,iw„ + ). (40) 

Equation (|38p now reduces to the quantum transport 
equation 

[iv n + e(k)/tl- e{k + )/h 

+IT (k, k + , iw„ , iio n+ )] $ a (k, k + , iw„ , ioj n+ ) 

-Za [£(k,iw n ) - g(k + ,iuj n+ )]p a (-q) 

= -iII 4 (k,k + ,ia;„,iw„ + )$ Q: (k,k + ,ia;„,iw„ + ). (41) 

Equation (|4"Tj) is the self-consistent equation for the auxil- 
iary electron-hole propagator, in contrast to the ordinary 
transport equation, equation (1471) , which is formulated in 
terms of the induced density 8n(k, q, iv n ) defined by 

Sn(k, q, iu n ) = - £ $ a (k, 

(42) 



The conditions under which (|4Tj) reduces to the ordinary 
transport equation are discussed below. 

To illustrate the role of the electron-EM field vertex 
renomalizations in (|4ip , it is useful to define the auxiliary 
single-electron propagator (?(k, ko n ) as the single-particle 
propagator which satisfies the following Dyson equation 

[ifej„ - e(k) + n - KE(k, iw»)]0(k, iu n ) = h. (43) 

The forward scattering processes cancel identically out 
in S(k, iw n ) as a result of the charge continuity equation. 
The fact that the auxiliary propagator (?(k, iw n ) has the 
quasi-particle pole even for relatively strong scattering 
processes is consistent with the fact that the residue asso- 
ciated with the quasi-particle pole in the exact propaga- 
tor Q(k,uu n ) is proportional to tb/t a and vanishes in the 
strictly ID limit if, -> 0. 30 This interesting observation is 
in the background of a variety of approximate treatments 
of various pseudogapped states. In these approaches, the 
forward scattering contributions to the FFCF .F(q', \v m ) 
are neglected from the outset. 

The expression for the intraband conductivity which 
correctly describes the u> « Drude limit, but which is 
not necessarily correct in the U) = Thomas-Fermi limit, 
is the following 

<r (q, w n ) = yJ2 ^(kK(-q)^: £ 

kcr iuj n 

x g(k,iu») - g(k+,iuj n+ ) 

\hu n + hU(k, k + , icj n , icj n+ ) + e(k) - e(k+) 

This is nothing but (|52"|) for the ideal intraband conduc- 
tivity in which irj is replaced by II(k, k + , iuj n , iw n+ ). 

B. Memory-function approximation 

In the weak coupling limit, the self-energy 
n(k, k + , i(jj n , vjJ n +) is a non-singular function of 
two fermion frequencies and in the simplest approxima- 
tion can be written in the form II(k, k_|_, iu> n , iu> n +) ~ 
M a (k, w n ) = M r a (k, \v n ) + iM* (k, xv n ) . M a (k, xv n ) 
is called the memory function^r— The momentum 
distribution function 

n W = -^J2 g ^n)e iu ^ (45) 

appearing in (|44[) in the memory-function approximation 
reduces in the weak coupling limit to the Fermi-Dirac 
distribution function f(k). 

We now define the electron-boson coupling function 
A Q (k, uj), the mass enhancement factor to* Q (k, uj)/m, the 
intraband relaxation rate r la (k, oj) and the collision in- 
tegral I(k, uj) by 

A a (k,w) = AfJ(k,w)/w, 

m aa( k > W ) = (1 + A "( k > w))m, 

r lQ (k, w ) = {m/m* aa (k,Lj))M l a (k,u), 

I(k,u>) = -T la (k,u)Sn(k,q,L)). (46) 
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Summation of (1411) over iw„, together with analytical con- 
tinuation \v n — > oj + vq and multiplication by E a (q,oj), 
gives the ordinary transport equation^^i^S 



-i[huj + E(k) -E(k + )]Sn(k, q,w) 



X^^E a (q,u;)=hI(k,uj) 



(47) 



where E(k) = e(k)/(l + A a (k, u = 0)). 

In solving (|47| it is convenient to follow the usual text- 
book procedure^ This equation can be solved either for 
J£(q, to) or for p c (q, uj), where 



(48) 



are, respectively, the induced charge and current densi- 
ties of conduction electrons, and the Sn.i(k, q, uj), i = 0, 1, 
are the contributions to 5n(k, q, uj) which are propor- 
tional to [i>°(k)]\ For further considerations of the 
random-phase approximation in section 8.1, the sec- 
ond route is more appropriate choice. In this case, 
we consider the thin slab with the longitudinal polar- 
ization of the electric fields, introduce the scalar fields 
V tot (q,uj) and V cxt (q,uj), by the relations V tot {q,uj) = 
(i/q a )E a (q,uj) and V cxt (q,uj) = (i/q a )E 0a (q,uj), 
and use the definition of the dielectric susceptibility 
X intra (q,w): p c (q,uj) = x intra (q, ^)^ tot (q, w ), «r 

[1/ X intra (q,c) - 47r/gV(q, W ) = V e *\q,u>). (49) 

It is readily seen that the intraband conductivity is of the 
form 



a^(q,uj) = ^X intra (q,^) = ie 2 ^E 

"at " i,„ 



V m* Q (k,w) 



[<(k)] 2 (-o/(k)/a £ (k)) 

> + iTi Q (k, uj) - q%v%Qs)/u ' 



(50) 



In the Drude ui rs limit, we set q a in the denominator 
equal to zero, with Ti a (k,uj) f=s (ri Q (k, w))fs — Fx a (u)) 
and m* aa (k,uj) (m* Q (k, w))fs = to obtain 



intra 



(q^) 



UJ + iTia(w) 



(51) 



an expression known as the generalized Drude formula. 
Here, (. . .)fs means the average over the Fermi sur- 
face. The expression (|5ip remains valid when more 
than one band intersects the Fermi level, provided that 
T^ L (ui) sa ri Q (w). The generalization to the case with 
the relaxation rates T^ l (uj) dependent on the band index 
L is straightforward. 

Although we have arrived at (l4"TT) and (|4T|) by consider- 
ing in the auxiliary self-energy S(k, iui n ) only the scatter- 
ing by boson modes, these expressions are of a more gen- 
eral validity, namely all relevant scattering mechanisms 



have to be included in II(k, k + , iui n , \ui n +)- Therefore, we 
can use (|4"Tj) to study the effects of local and short-range 
electron-electron interactions on electrodynamic proper- 
ties of strongly-correlated low-dimensional systems or to 
study similar effects associated with the scattering by 
different periodic CDW/SDW potentials. 



VII. INTERBAND CONDUCTIVITY 

In analogy to (|44p . we can write the bare interband 
conductivity in the form 



<T'°^n) = ~ E E^E^ L '( k )^ L ' L ( k ) 

Qb{k, \u>n) - GL'{k+,iuJ n+ ) 

ifw n + hU LL ,(k, k + ,iw„,iw n+ ) - e L ,{k + ) + e L (k) ' 

(52) 

with the electron-EM field vertex corrections 
neglected, resulting in II(k, k + , iu> n , ioj n +) ~ 
E(k, iui n ) — £(k+, kj„+). The electron- hole self-energy 
II LL /(k, k + ,kj„,iw n+ ) is replaced by M^ L ' (k,iv n ), 
in the memory-function approximation, and by 
in the relaxation-time approxi- 



ir^ (k) 



ir 



LL' 



mation. 

In those multiband systems in which the interband 
dipole vertex P^ L (k) is purely imaginary, the lo- 
cal dipole moment vanishes and there are no local 
field effects.— The interband conductivity is given by 
a™^ er (q, w) = <7^* er '°(q,u;). This result can be equally 
well understood to represent the solution to the RPA 
equation 

[1/ X (q, u) - ^/q 2 ]p(q, uj) = V cxt (q, w) (53) 
for the total dielectric susceptibility xiQi 1 ^) — 

E i= i„tra,interX l (q,w), where x'fa") 

(9«/iw)^ Q (q- w )- 

In the opposite case, the interband dipole vertex is 
real and its dispersionless part is directly related to the 
intracell dipole transitions. Particularly interesting are 
the examples in which P^ L (k+ , k) s» P^ L is completely 
local. In this case we can use the procedure from the 
following section to determine the effects of the Lorentz 
local fields. 

At this point caution is in order regarding the 
relaxation-time approximation. It is well known that in 
those systems in which the damping energy T^ L (k) is 
comparable to the interband threshold energy (the pseu- 
dogapped state of the ordinary CDW systems is an ex- 
ample) the relaxation-time approximation of (1521) overes- 
timates the in-gap part of the interband conductivity. In 
this regime, 0-™„ er '°(q, uj) is found to be better described 
by the approximate textbook expression 38,39 



P Q "'(k)| 2 [/ L (k)-/ L ,(k + )] 



V ^ foj + ihT 2a (k) -e L >(k+) + e L (k) • 



(54) 
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FIG. 9: The mean-field single-particle conductivity of the 
Peierls CDW model along the highly conducting direction x, 
at temperature T — 160 K. The bare electron dispersion is 
e(k) = — Yl a =x y % a c °sk a a a , with 2t a = 0.5 eV and 2tt = 
50 meV. The mean-field transition temperature is assumed 
to be T£ = 200 K and A(160K) = 21 meV is the CDW 
order parameter at temperature T = 160 K. The damping 
energies are 5KTi = KC 2 = 30 meV at T < T$, and hT'l = 30 
meV at T > Tp. The dot-dot-dashed line shows the inverse 
square-root singularity in a'^ CI ,sp (oj) associated with the limit 

hv 2 -»• o. 



In order to briefly illustrate our results, the typi- 
cal mean-field single-particle conductivity of the Peierls 
CDW model at temperatures well below the mean-field 
transition temperature Tp is shown in figure 9. The 
dashed and the dot-dashed lines represent, respectively, 
the interband contributions calculated by using (1541) and 
the relaxation-time version of (|52p. Notice the differ- 
ence between the two curves at w « 0. The solid line 
is the total single-particle conductivity which is the sum 
of the contribution of thermally activated charge carriers 
in (|5ip and the interband contribution (|54|). This result 
agrees reasonably well with the experimental spectr 
measured at temperatures well below the CDW transi- 
tion temperature. The detailed discussion of the mean- 
field approximation as well as of the effects of fluctuations 
in the order parameter is given in Article II. 

In the next section, we shall examine the Lorentz local 
field corrections in the relaxation-time approximation, 
starting with <7^ er '° '(q, uS) given by (|54p. The damp- 
ing energy HT 2a (k) rs fiT 2a is assumed to be small in 
comparison with the interband threshold energy. 



VIII. LOCAL FIELD EFFECTS 

Most Q1D systems are characterized by a large number 
of atoms in the unit cell (for example, in the blue bronze 
K0.3M0O3 we have 86 atoms per unit cell). The bands 
in the vicinity of the Fermi level are thus built of differ- 
ent molecular orbitals of the well-defined symmetry with 



the intermolecular distance large when compared to the 
first-neighbor atom-atom distance. Therefore, the bands 
are narrow and the typical interband threshold energy is 
in the infrared part of the spectrum. In most cases of in- 
terest, the low-energy interband electron- hole excitations 
are local and correspond to dipolar transitions between 
one molecular orbital of the s symmetry and another of 
the p a symmetry. Consequently, to determine the proper 
structure of the total conductivity tensor er aQ (q, w) in 
this case, it is important to take into account the local 
field effects too. 

This is not the only Q1D case in which the local field 
corrections are important. In the ordered state of the 
strong bond-CDW systems the interband dipole excita- 
tions are also local. Finally, the ordinary sp a models 
with the interband threshold energy in the visible part 
of the spectrum also have tr™* Gr '°(q, to) of the form (|54p . 



with the dipole vertex P£' L (k) ss P£' L real 



25.26,39 



A. Multicomponent random-phase approximation 

To demonstrate the usual procedure for dealing with 
the local field corrections j2£r— we consider the model 
with the intraband currents and with one type of bound 
charges. We assume that the contribution of the bound 
charges is related to the intra-atomic or intramolecular 
dipolar transitions in erj^ er >°(q, lj). In section 8.2, we 
shall briefly discuss the case of infrared-active optical 
phonons, as well. In this section, we make a simple mod- 
ification of the notation. For example, the label intra is 
replaced by c and the label inter by b (for conduction 
and bound charges). 

The bare intraband dielectric susceptibility is assumed 
to be of the form 



X c '°(q, W ) ee Xii(q,^) = (&Mo*£(q,uj), (55) 



— ,,0 



c,Q/ 



and the bare interband 



with cr^°(q,w) given by 1-5 
dielectric susceptibility is 



X b <°(q,^) = 5&(q,a0 = (ql/i^h^u), (56) 

with o-k'°(q,w) given by (|54). 

The electrons are assumed to be subjected to an ex- 
ternal longitudinal electric field E (r, t) represented by 



the scalar potential F ext (r, t). pi(q, 



p c (q, uj) and 



p 2 {(\, w) = — iq a P b (q, lu) are the Fourier transforms of 
two types of macroscopic charge densities. /0 c (q, ui) de- 
scribes the monopole-charge density of conduction elec- 
trons and P^(q, w) the dipole-charge density of bound 
charges. The two-component random-phase approxima- 
tion describes the long-range screening in such a case. 
According to figure 10, the RPA equations can be shown 
in the form 

[1/X?i(cbw) - V u ]pi(q,w) ^ V 12P2 (q,u) = V cxt (q,uj), 

+ [i/x2 2 (q,w) - v 22 ] P2 (c l ,u) = r= xt (q,^)- 

(57) 
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FIG. 10: The RPA equations for the screened susceptibilities 
Xii(q,w) (a) and X22(q,oj) (b). 

Here, the Xui^ 10 ) stand for two bare dielectric suscep- 
tibilities in ((SSJ) and (|55]). and Vij = Vy(q) = -Nj/q 2 
are the Fourier transforms of different long-range inter- 
actions. The numbers Nu = —4ir, N12 = N21 = — 47r 
and N22 = — 47r + Att/3 correspond to the monopole- 
monopole, monopole-dipole and dipolc-dipolc interac- 
tions, respectively^ Nu — N 12 = N21 = — 47r are easily 
shown to represent the common depolarization factors 
for the longitudinal q ~ fields^ It is not hard to verify 
that, after neglecting the Lorentz term 47r/3 in N22 (i.e., 
when N22 = Nu), equations (|57)) reduce to (|53)) . 

Equations (|57|l may be rewritten now in the alterna- 
tive form, in terms of two dipole densities Pj a (q, ui) = 
(V?«)Pi(q,w), to give 

(l/a° u )P la = E 0a + N n P la + N 12 P 2a , 

(l/a° 2 + N n - N 2 2)P2a = E 0a + N 21 Pi a + N 11 P 2a , 

(58) 

with the implicit dependence on q and ui. Here 

e&(q J w) = (-l/«£)x&(q,w) (59) 

is the common notation for two bare polarizabilities. It 
must be emphasized that the expressions on the right- 
hand side of (|58p represent the macroscopic electric field 
E(q, ui), equation (|3T))) . which is the sum of the exter- 
nal electric field and two depolarization fields originat- 
ing from the conduction electrons and from the induced 
dipole moments of bound charges. 14 This is a useful ob- 
servation which helps us to rewrite ([55} in the final form 

[l/a c aa (q,oj)]j c a (q,^) = E a (q,u), (60) 

[lK Q (q, W )]P Q b (q,w) = E a (q,w), (61) 

with 

E a (q,u) =^a(q,w)-47r[Pi (q J w) + i^„(q,w)]. (62) 

For the thin slab with the longitudinal polarization of 
the fields, the field Eo(r,t) is identical to the dielectric 



displacement field D(r, t) and the relation (|62|) leads to 
the usual form of the longitudinal macroscopic dielectric 
function. If other high-energy on-site dipolar excitations 
are also included, the total dielectric function becomes 

£QQ (q,w) = £™(q,w) + ^a*^) +4^ a (q,c). 

(63) 

Here, a u (q,w) = a a n (q,uj) = (i/u;)o£ a (q,u;) and 

a22(q '" } = l-(47r/3K 2 (q,a,) = (64) 

are two renormalized polarizabilities. 

In the notation of the previous sections, this result can 
be expressed by (fT0|) and (fTTj) . We can therefore conclude 
that the contribution of the monopole charges to the con- 
ductivity tensor is decoupled from the contribution of 
the dipole charges. The latter contribution is affected by 
the Lorentz field in the usual way, giving also the possi- 
bility for the collective contribution in c™^ er (q, ui) when 
1 + (47r/3w)Im{<^°(q, u)} = 0. 

B. Infrared-active optical phonons 

It was implied by our discussion of the quantum trans- 
port equation that ([SJ, together with the density opera- 
tors <j3j) and (|23|). is able to capture all relevant aspects 
of the electron-phonon coupling in the dynamical con- 
ductivity tensor of multiband electronic systems, at least 
for the scattering by acoustic and Raman-active optical 
phonons. In order to treat carefully infrared-active op- 
tical phonons, we must also include the direct coupling 
between phonons and EM fields in H cxt and take into 
account the corresponding local field corrections. 

The generalization of (|63| to such a case with multiple 
dipolar bound-charge contributions is very simple. For 
our present purposes it is simplest to replace cvS^q, w) by 
the sum of all bare dipolar bound-charge polarizabilities. 
In the two-band case with one infrared-active phonon 
branch, we obtain 

ag 3 (q,w) = (-l/^)x intor,0 (q,w) +a° h (q,w), (65) 

with a£ h (q,w) = (e 2 /A^V u )[w2 q - w (a; + i 7 )]- 1 being the 
bare phonon polarizability. 

IX. CONCLUSIONS 

In this article, we have derived the dynamical con- 
ductivity tensor in a general multiband electronic model 
with sizeable boson-mediated electron-electron interac- 
tions. Starting with the gauge-invariant form of the cou- 
pling between valence electrons and external EM fields, 
we have derived the general form of the self-consistent 
Bethe-Salpeter equations for the auxiliary intraband and 
interband electron-hole propagators. It is demonstrated 
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that the intraband Bethe-Salpeter equation is identical 
to the usual quantum transport equation and that it re- 
duces to the ordinary transport equation in the weak 
coupling limit. In this way, the present response the- 
ory of multiband electronic systems is connected to the 
generalized Drude formula for the intraband conductiv- 
ity widely used in experimental analyses as well as to 
different Fermi liquid theories based on the Landau-Silin 
transport equations. In order to illustrate the formal- 
ism presented here, we shall give in the accompanying 
article^ the detailed description of the CDW mean-field 



approximation and the effects of fluctuations in the order 
parameter on the dynamical conductivity of the Peierls 
CDW systems. 
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